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LAN PROPERTY FOR STOCHASTIC DIFFERENTIAL EQUATIONS 
WITH ADDITIVE FRACTIONAL NOISE 
AND CONTINUOUS TIME OBSERVATION 

EULALIA NUALART AND SAMY TINDEL 


Abstract. We consider a stochastic differential equation with additive fractional noise of 
Hurst parameter H > 1/2, and a non-linear drift depending on an unknown parameter. We 
show the Local Asymptotic Normality property (LAN) of this parametric model with rate 
\/t as r —>■ oo, when the solution is observed continuously on the time interval [0,r]. The 
proof uses ergodic properties of the equation and a Poincare type inequality. 


1. Introduction 

Let R be a </-dimensional fractional Brownian motion with Hurst parameter H > 1/2. 
Let us recall that R is a centered Gaussian process defined on a complete probability space 
(G, R, P). The law of R is thus characterized by its covariance function, which is defined by 

R,;i^E[R/R/] + l{o}(*-j), s,teR. ( 1 ) 

The variance of the increments of R is then given by 

E[|Ri-R:r] = * = ( 2 ) 

and this implies that almost surely the fBm paths are y-Holder continuous for any 7 < R. 
In this article, we consider the following RGyalued stochastic differential equation driven 

byR; 

yt = yo+ b(Ys-,9)ds+ te[0,r]. (3) 

Jo ,=i 

Here ijq G is a given initial condition, R = (R^,..., R'^) is the aformentioned fractional 
Brownian motion (fBm), the unknown parameter 6 lies in a certain set 0 which will be 
specified later on, { 6 (-; 6), 6 & 0} is a known family of drift coefficients with 6 (-; 0 ) : ^ 

and (Ti,..., (Trf G are assumed to be known diffusion coefficients. 

There has been a wide interest in drift estimation for stochastic equations driven by 
fractional Brownian motion in the recent past, partly motivated by inverse problems in a 
biomedical context [20]. However, notice that the existing literature on the topic mainly 
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focuses on the case where the dependence 6 i-A- b{x] 6) is linear and all coefficients are real¬ 
valued. In this situation least squares and maximum likelihood estimators (MLE) for the 
unknown parameter 6 can be computed explicitly, and numerous results are available: the 
case of continuous observations of Y and where the drift coefficient is linear in both 6 and x 
(that is, fractional Ornstein-Uhlenbeck type process) has been studied e.g. in [2, 16, 18, 23], 
either in the ergodic [9 < 0) and non-ergodic {9 > 0) case. See also the monograph [30] 
and the references therein. Results on parameter estimation based on discrete observations 
of Y in the linear case can be found e.g. in [4, 35]. The case of a dependence of the form 
{9^x) i-A- 9h{x) is handled e.g. in [21, 33]. 

However, the case of a general multi-dimensional coefficient b{9, x) in equation (3) is also 
a very natural situation to consider, though we are only aware of the articles [6, 27] giving 
some positive answers for such a dependence. The current contribution has thus to be 
thought of as a step in this direction. Indeed, our main aim is to show that the model given 
by equation (3) satishes the Local Asymptotic Normality property (LAN) with rate ^/T as 
r —)■ oo, when the process Y is observed continuously on the time interval [0, r]. 

Before we proceed to a specific statement of our results, let us describe the assumptions 
we shall work with, which are similar to the ones used in [27]. We start with a standard 
hypothesis on the parameter set 0: 

Hypothesis 1.1. The set 0 is compactly embedded in Md for a given q > 1. 

In order to describe the assumptions on our coefficients b, we will use the following notation 
for partial derivatives: 

Notation 1.2. Let / : x 0 —)■ M 6e a function for pi,p 2 > 0. Then for any 

tuple (A,... ip) G {1,..., d}P, we set for Analogously, we use the notation 

f ; whcrc (ii,... ip) is a tuple in {1,..., q}P. Moreover, we will write dxf 

resp. dof for the Jacobi-matrices {dx^f ,..., dx^f) and {de^f ,..., do^f). Finally, for the sake 
of simplicity, we denote by (•, •) the Euclidean scalar product in or Mf", and by \ ■ \ the 
corresponding Euclidean norm. 

Let us now state the linear growth plus inward assumptions on our drift coefficient ensuring 
ergodic properties for the process Y : 

Hypothesis 1.3. (i) There exists ci > 0 such that for every x,y & and 9 E Q we have 

{b{x; 9) - b{y; 9), x-y) <-a\x- y\‘^. 

(ii) ITe have b G x 0;M'^), with dxb and d‘lJ) uniformly bounded in {x,9). 

(hi) We have b := dgb G C^’°(M'^ x 0; x M^), with dxb and d^J) uniformly bounded in (x, 9). 

(iv) There exists c > 0 such that for every x,y eMA and 9, 9 q E 0, the following Lipschitz 
and growth conditions hold: 

\b{x-,9)\ < c(l |x|), \b{^;9)\ < c(l |x|) 

\b{x-,9) -b{x;9o)\ < c\9 - 9o\{l + |x|), \b{x;9) -b{y;9)\ < c\x - y\. 

Furthermore, we suppose that the coefficient a fulhlls an invertibility condition of the 
following form: 
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Hypothesis 1.4. The number d of driving fBm equals the dimension of the state space 
for Y. In addition, denoting by a the d x d matrix with columns ai,, ad, we assume that 
a is invertible. 


Notice that Hypothesis 1.3 and 1.4 entail the following result (see next section for more 
details); for a given 6 * G 0 the solution to equation (3) converges for t —)■ cxo almost surely 
to a stationary and ergodic stochastic process {Yt,t > 0 ). 

Let us now introduce the concept of LAN property in our context. Towards this aim, for 
any 6 E Q and X < H, we denote by (resp. PJ) the probability laws of the solution to 
equation (3) in the spaces (resp. C^([0, r]; Moreover, we assume that the 

process Y is observed continuously in [0,r]. In this context, the dehnition of LAN property 
takes the following form; 

Definition 1.5. We say that the parametric statistical model {Pe, d G 0} satisfies the LAN 
property at 6 E Q if there exist: 

{i) A q X q invertible matrix 

(ii) A q X q positive definite matrix r(6'), 

such that for any u E the following limit holds true as r ^ oo: 

log «W(0,r(9)) - ivr(9)«, 

where J\f{ 0 ,r{ 6 )) is a centered q-dimensional Gaussian random variable with covariance 
matrix r(0). 


The LAN property is a fundamental concept in asymptotic theory of statistics, which 
was developed by Le Cam [24] . The essence of LAN is that the local log likelihood ratio is 
asymptotically normally distributed, with a locally constant covariance matrix and a mean 
equal to minus one half the variance. The main application of the LAN property is the 
following Minimax Theorem. It asserts (if LAN holds true) an asymptotic (minimax) lower 
bound for the risk with respect to a loss function, for any sequence 9^ of estimators of 9. 


Theorem 1.6 (Minimax Theorem). [17, Theorem 12.1] Suppose that the family of proba¬ 
bility measures {P 0 )e^e satisfies the LAN property at a point 9. Let (9 t-)t>o be a family of 
estimators of the parameter 9 and / ; —)■ [0, +cxo) be a loss function such that 1(0) = 0, 
l{x) = l{\x\) and l{\x\) < l{\y\) if |x| < \y\. Then we have: 


lim lim inf sup E^/ 

(5-s-o T^oo 




where C{Z) = (9)-^). 


> E. [/ (Z)], 


( 4 ) 


Observe that when the loss function I in Theorem 1.6 is the quadratic function, the lower 
bound in (4) is simply the trace of r(0)“^. A sequence of estimators that attains this 
asymptotic bound is called asymptotic minimax efficient. Therefore, Theorem 1.6 opens the 
way to a theory which mimics the concept of efficient estimator from the Cramer-Rao lower 
bound. 
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As one can see, the LAN property is an important tool in order to quantify the iden- 
tifiability of a system. The aim of this paper is thus to show the following result for our 
equation (3). 


Theorem 1.7. Assume Hypothesis 1.3, and recall that stands for the probability law of 
the solution to equation (3) in the space C^([0, r]; M'^). Then, for any 0 G 0 and u & 
fixed, as T ^ oo, we have: 

dPZ,_^ £/p \ 1 

log -> «’'A/'(0, r(«)) - (5) 

6 

where the matrix r( 6 ') is defined by: 


r( 0 ) := 


Eg[{b{Yo;e)-b{Yrp,e)y 


cr 


-DT^-l 


fb{Yo;e)-b{Yr,;e))] 


1/2+H 1/2+H 


dridr2, ( 6 ) 


In (6), recall that Y is the ergodic limit of our process Y and that we have set b for the 
Jacobian matrix dgb. 


Though our result encompasses a wide range of coefficients b, it is worth illustrating it 
on the simple linear case, that is for the real-valued fractional Ornstein-Uhlenbeck process 
corresponding to b{x] 6) = —9x and a = 1. In this case we have: 

Yt = —6 [ Yg ds + Bt, f G [0, r], 6 > 0. 

Jo 

Let 6r be the MLE of 6 from the observation of Y in [0,r]. Then, it is well-known (see 
[3, 18]) that 

lim 9^ = 9 a.s. 

r—>-OD 

and that 

/:(Pe) - lim y/f{9r -9)= M{0, 29). (7) 

r^oo 

In particular, this last result already suggests that the rate of convergence for the LAN 
property in this case is as mentioned in [30, p.l62]. Our result confirms this intuition, 

and one can see that the MLE reaches this optimal rate of order However, observe 

that the identification of ( 6 ) and (7) is far from being clear, leaving open the question of 
efficiency of the ML estimator in the fractional Ornstein-Uhlenbeck case. 

The LAN property for stochastic processes has been widely explored in the literature. To 
mention a few references close to our contribution, the case of ergodic diffusion processes 
observed continuously (with an unknown parameter on the drift coefficient) is dealt with in 
e.g. [22, Proposition 2.2]. The proof follows from Girsanov’s theorem, ergodic properties 
and the central limit theorem for Brownian martingales. The rate of convergence achieved 
in [22, Proposition 2.2] is ^/T like in our Theorem 1.7, with a function P(0) given by: 

P(0) = Eglb(Y;9)^a-\V)^a-yY)b(Y;9)], 

which should be compared to our expression ( 6 ). The case where the Brownian diffusion 
process is observed discretely was solved in [ 11 ] using an integration by parts formula taken 
from the Malliavin calculus. Let us also mention that the LAN property is obtained in [7] for 
some discretely observed fractional noises including fractional Brownian motion. This last 
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result is achieved by a direct expansion of log-likelihood functions, plus a thorough analysis 
based on properties of Toeplitz matrices. 

Finally let us say a few words about the strategy we have followed in order to handle 
the case of equation (3). In order to compare likelihood functions for different values of 
the parameter, we take the obvious choice of applying Girsanov’s theorem for fractional 
Brownian motion (following the steps of [26]). We are then left with two technical difficulties 
in order to get asymptotic results: (i) Handle the singularities popping out of the fractional 
derivatives in the Girsanov exponent, (ii) Get ergodic results in Holder type norms for our 
process Y. Let us also mention that GLTs for martingales, which were an essential tool in 
the diffusion case, are unavailable in our fBm situation. We sort out this problem by means 
of Malliavin calculus techniques. Specihcally, we shall derive some concentration properties 
for the Girsanov exponent by means of a Poincare type inequality (see relation (16) below). 
This means that some Malliavin derivatives will have to be conveniently upper bounded, 
which is a non negligible part of our efforts. 

Our paper is structured as follows: Section 2 is devoted to preliminary results, including 
Malliavin calculus background and the derivation of most of the concentration properties 
for the process Y needed later on. Then we prove Theorem 1.7 in Section 3, following the 
strategy described above. 

Let us close this introduction by giving a set of notations which will prevail until the end 
of the paper. 

Notation 1.8. We use the following conventions: for 2 quantities a and b, we write a 
if there exists a universal constant c such that a < cb. In the same way, we write a ^ b 
whenever a <b and b Y a. If f is a vector-valued function defined on an interval [0, r] and 
s,t ^ [0,t], 6fst designates the increment ft — fg. 

2. Auxiliary Results 

In this section we hrst recall some basic facts about stochastic analysis for a fBm B, and 
also about Young integrals. Then we shall derive some pathwise and probabilistic estimates 
for equation (3) under the coercive Hypothesis 1.3. 

2.1. Stochastic analysis related to B. The ergodic properties of equation (3) are ac¬ 
curately encoded by the fBm representation given in [15], which goes back to the original 
paper [25]. We present this construction here for a one-dimensional fBm, for sake of simplic¬ 
ity. Obvious generalizations to the d-dimensional case are left to the reader. 

2.1.1. Wiener space related to B. Let IT be a two-sided Brownian motion, and H a Hurst 
parameter in (0,1). We consider the process B dehned for t G M by: 

Bt = CH [ [dWt+r - dWr] (8) 

Jr- 

= Ch\ [ [(-(r - f))^-^/2 - (-r)^-i/2j - [ {-{r-t))^-^/^dWr 

iJ-oo Jo 

Then B is well-defined as a fBm, that is a centered Gaussian process with covariance given 

by (1). 
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The process B is closely related to the following fractional derivatives: for a G (0,1) and 
(p G we set 

(p] ^ = Ca / and = Ca / (pt-rr°~^dr. 

With this notation in mind, the following proposition identifies a convenient operator which 
transform W into B as in (8). 

Proposition 2.1. For w G C^(IR) and H G (0,1), set 

[KHW]^ = ch [wt+r - Wr] dv. 

Jr. 


Then the following holds true: 

(i) There exists a constant ch such that 

for H > 
for H < 

(ii) For H > 1/2, Kh can he extended as an isometry from L^(M) to 
(hi) There exists a constant ch such that we have Kfj^ = chKi_h- 


[Khw\ = 


-Ch ([I+ - [1+ , 

-Ch [[D]!'^~^w]t - 


1 

2 

1 

2 ‘ 

1/2 


(L2(M)). 


The natnral Wiener space related to our fractional Brownian motion indexed by M, as 
represented in (8), is given as follows (see [15] for further details): one can take B as the 
underlying space, where 


B 


f e C(M; M"); 


\ft 


1 + \t\ 


< oo 


For this choice of B, equipped with the law P of our fractional Brownian motion B, there 
are various ways to express the generic form of an element of the Cameron-Martin space 
Fi. Among those we will choose the following one: h is an element of Fi iff there exists an 
element Xh in the first chaos such that: 


ht = F.[BtXhi and ||h||^ = (9) 

Moreover, one can check that {B,Fi,'P) is an abstract Wiener space. Notice that the space 
Fi has to be distinguished from the space Fi usually considered (e.g in [28]) for the Malliavin 
calculus with respect to B. For the Wiener process IF, we would have for instance Fi = 
L^(M), while Fi has to be identified with the Sobolev space Fr^(M). 

Let us recall a definition concerning variation norms which can be found e.g in [10]. For 
0 < s < f < r and a function / : [s, t] —)■ M'^, we set: 

||/i||i-var;p,f] = sup ^ fuu+i \ ] kl,,* partition of [s,t] 

i. ti Gllst 

With this definition in hand, we can label the following property of elements of Fi for further 
use: 

||/l||l-var;M < (^ “ s)^ ]| h]| ^. (10) 
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Next we give a representation of the Canieron-Martin space on a finite interval [0,r], 
namely ?^([0,r]), which can be found in [ 8 ]. To this aim, let us recall the dehnition of 
fractional integrals of order a G (0,1] on a hnite interval: 

= f\t- nr-'.Udu, and [/“_/], = ^ /"(a - ()“’'/.<;«■ (H) 

r(a) Jo r(a) Jt 

With this additional definition in hand and assuming H > 1/2, a function h lies into the 
space ?/([ 0 , r]) if it can be written as 

h = krf, with krf = and /eL^([ 0 ,r]), ( 12 ) 

where we have set k = H — 1/2. Related to this representation, we have ||h||-^Qo,T]) = 

||/||l 2 ([o,t]) whenever h = K^f. Notice that this space will play a prominent role in our 
estimates on finite intervals, associated to another space called ?f([ 0 ,r]) which pops out in 
Malliavin calculus computations (as mentioned above for Ji and ?/). Indeed, if h G ”^([0, r]) 
admits the representation ( 12 ), then for a smooth function g on [ 0 , r] we have: 

f gsdhs= [ gss'^[lQ^{x~'^f)]sds = f [K*g]sfsds, (13) 

Jo Jo Jo 

where we have set: 

K;g = x-^U4x^g), 

and we recall that Ir- is defined by ( 11 ). and where the integral gg dhg is understood in the 
Young sense. The closure of smooth functions under the inner product {K*gi, R"/fi' 2 )L 2 ([o,r]) 
is the space which is called 'H([0,r]). We shall resort to the following formulae (for which 
we refer e.g to [ 8 ]) concerning 'H([ 0 ,r]) in the sequel: 

\\9\\n([o,T]) = CH 9s^\si-S2f^~^gs2dsids2, ||^|k([o,r]) < (14) 

J[0,r]2 

The following relation between norms in "H and will also be needed in the sequel. 

We include its proof here for sake of completeness, because we haven’t been able to track it 
down in the literature. 

Proposition 2.2. Let H > 1/2 and r > 0. Consider the space 'H([0,r]) defined by (12) and 
the usual Sobolev space kT^’^([0, r]). Then for all h G 'H([0,r]) we have: 

\\h\\wr2([o,r]) < ch \\h\\ii. (15) 


Proof. Start from the representation (12) for a typical element h G K([0, r]), and recall that 
we have set k = H — 1/2. Then we have: 

J-s hg ds , II-^11 l 2 ([o,t]) 1 

Consider now a generic element i G L^([0,r]) such that ||^||l2([o,t]) = 1- We resort to (13) in 
order to get: 


WL 2 ([ 0 ,t]) = sup 




< ii^;^ii 


L2([0,r]) 


L2([0,r]) 


\\K*J\\ 


L2([0,r]) 


mO,T]) ’ 
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where we recall that / comes from representation (12). Furthermore, invoking our bound (14), 
we end up with: 



which proves our claim. 




n{[ 0 ,r]) 


\\h\\n([0,r]) ^ 


□ 


2.1.2. Malliavin calculus tools. We refer to [28] and [34] for basic definitions concerning 
Malliavin derivatives D and Sobolev spaces related to B. Our concentration arguments 
for the solution to (3) will heavily rely on the following Poincare type inequality borrowed 
from [34, page 76]. 

Proposition 2.3. Let F : B ^ M. be a functional in D^’^. Then, 

E||F-E|F]|]<|E[||flF||fl], (16) 

where we recall that D denotes the Malliavin derivative. 


Next we prepare the ground for our Girsanov transform by relating a general drift with 
respect to a fractional Brownian motion with its counterpart with respect to W (we con¬ 
tinue in the one-dimensional setting for the sake of simplicity). This proposition is a slight 
extension of [15, Lemma 4.2]. 


Proposition 2.4. Let H G (1/2,1), he two paths in and suppose that 

F = Khw^ for j = 1,2, where Kh is defined at Proposition 2.1. We also assume that for 
t>t],F and 6 ^ are linked by the relation: 


bl = bl, and b'f = b] + / gb,s ds, 


for a function in L‘ 


Then we also have: 
wl =wl+ / g^^s ds, with 51 ^ = D^~^^‘^gb. 


Proof. By [15, Lemma 4.2, relation (4.11a)] and [29, relation (4)], 
d /■* 


_ 9b,u 

9w,t (t - m)^^-V2 

which finishes our proof. 


du = ch 




T^H-I/2 

D+ 9b 


□ 


The previous proposition yields a version of Girsanov’s theorem in a fractional Brownian 
context, inspired by [26] (see also [29]); 

Proposition 2.5. Let B he a fractional Brownian motion and g he an adapted process in 
/;J^~i/^(L^(M)). We set F) = j^ggds and Qt = Bt + Vt for t > 0. Then Q is a fractional 
Brownian motion on [0,r], under the probability P defined by: 


dP 

dP |[o,T 
W-1/2 


= e 


-L 


where L = 


L= rp 

Jo 


H-1/2 


9l„ciir„ + i / [D"-''^g]ldu, 


(17) 


provided Dj^ ' g satisfies Novikov’s conditions on [0,r]. 
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2.2. Generalized Riemann-Stieltjes Integrals. We are focusing here on the case of a 
Hurst parameter H > 1/2 for sake of simplicity, so that all the stochastic integrals with 
respect to B should be understood in the Young (or Riemann-Stieltjes) sense. In order to 
recall the definition of this integral, we set 


||/||oo;M = sup |/(t)|, 
tG[a, 6 ] 


|/|A;[a,fe] 


sup 

s,t£[a,b] 


\t - s|^ 


where / : M —)■ M” and A G (0,1). We also recall the classical definition of Holder-continuous 
functions; 


C"(K6]iR”) = {/:|a,6|^ 


3;|a,6| + |/U;[a,b| ^ Oo| . 


Now, let / G C^{[a, b]] M) and g G C^{[a^ &]; M) with A -|- p > 1. Then it is well known that 
the Riemann-Stieltjes integral fs dgs exists, and can be expressed as a limit of Riemann 
sums. 

In the sequel, the only property on Young’s integral we shall use is the classical chain rule 
for changes of variables. Indeed, let / G G^([a, 6];M) with A > 1/2 and F G Then 

we have: 

FUy) - F{fa) = r F'if,) dfs, y G [a, b]. (18) 

J a 


2.3. Basic properties of solutions to SDEs. We first state the existence and uniqueness 
result for equation (3) borrowed from [15, Lemma 3.9]. 

Proposition 2.6. Under Hypothesis (1-3), there exists a unique continuous pathwise solution 
to equation (3) on any arbitrary interval [0,r]. Moreover the map Y : {yQ,B) G M'’* x 
C([0, rj; M'’^) —)■ C([0, rj; R'’*) is locally Lipschitz continuous. 

In addition, exploiting the integrability properties of the stationary fractional Ornstein- 
Uhlenbeck process, the following uniform estimates hold true. They are shown with more 
details in [12], and are obtained by comparing Y with the solution to the equation dUt = 
—Ut dt + dBt (much in the same way as in the next section 2.4). 

Proposition 2.7. Assume Hypothesis 1.3 holds true and let Y be the solution to equation (3). 
Then for any 0 G 0 and p>l there exist constants Cp, kp > 0 such that 

E [\Yt\P] < Cp, and E [|<5W/p] < kp\t - s\p^, 

for all s,t > 0, where we recall that we have set SYgt = Yt — Yg as mentioned in Notation 1.8. 

Let us now state a pathwise bound on the increments of Y, where we use a similar idea 
as in [19], in order to translate the moment estimates into pathwise bounds. 

Proposition 2.8. Assume Hypothesis 1.3 holds true and let Y be the solution to equation (3). 
Then for alls G (0, H) there exists a random variable G r\p>iLP{Q) such that almost surely 
we have: 

\Yt\<Z,{l+tf^, and \6Ygt\<Z,{l+tf^\t-sf-^ 
uniformly for D < s <t. 


( 19 ) 
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Proof. We focus on the bound for SYgt, the estimate on |h^| being shown in the same manner. 
Set 'j = H — e and consider n > 1. Thanks to Garsia’s lemma (see e.g. [13]) and recalling 
our Notation 1.8, for all e [0,n] we have where p > 1 and Af^y^p^n 

is defined by: 

r |<5w 


Af'y^p^n — 


'0 JO 


- vYp+‘^ 


\v — u 


p \p 

dudv j 


Furthermore, owing to Proposition 2.7 we immediately get: 

E'''” :£ 

under the constraint p > Set now — supngj.^ Then for all q > we have: 


E [ZJ] = E 


sup —^ 

nGN 


^ ^ E [|A/1 


<ji _ _ 1 

I2vhi<y^<oo, 


n=l 


^2eq 


n=l 


( 20 ) 


namely we have found G np>£-iL^(r2). Obviously this also yields G np>iL^(r2). Finally, 
replace n > 1 by [f] + 2 to obtain our pathwise bound (19). □ 

Let h be a generic element of PL (this space being defined by (9)), and set for the 
solution to our equation (3) driven by i? + h. We shall need an extension of the previous 
proposition to Y^, whose proof is left to the reader since it follows along the same lines as 
Proposition 2.8. 

Proposition 2.9. For h G PL and 9 E Q, let Y^ be the solution to: 

rt d, 


L* = v«+ i KP; i))ds + Y, {B’l +l‘i), te [o, t] . 

i=i 


( 21 ) 


>0 


Under the same assumptions as for Proposition for all e G (0, H) there exists a random 
variable Z^ G np>iL^(G) such that almost surely we have: 


\5Yf,\<Z,{l+tt{l + \\hU)\t-s 


\H-e 


( 22 ) 


uniformly for 0 < s <t. 

We next show that our Girsanov type theorem 2.5 can be applied to the solution to (3). 

Proposition 2.10. Assume Hypothesis 1.3 and let Y be the solution to (3), for a given 
6 E Q. Consider the d-dimensional process Qt given by: 

Qt= [ a-^b{Ys-,e)ds + Bt. 


Then Q is a d-dimensional fractional Brownian motion under the probability P^ defined by 
dPe — ^-L 
dPe 


l[o,- 


= e ^, with 


L= {a-^[Df-^^hiY-e)UdWu) + f- \a-^[Df-^^\Y-e)U^du. 


(23) 
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Proof. We first need to show that ; 6) is well dehned on [0, r], which follows since 

the fnnction u i—)■ b{Yu] 9) is {H — e)-Holder, and thns {H — l/2)-Holder. 

Second, we need to verify Novikov’s condition (see [9, Theorem 1.1]), that is, that there 
exists A > 0 snch that 


sup Eg 

4G[0,t] 


exp 


A 


cr 


-ir 


[D^-^^\Y-9)]fds 


< oo. 


(24) 


Towards this aim, we extend the dehnition of Y to M_ by setting h) 


O' 


-1 


[D^-^^\Y)l\ = CH 


a 


-1 


{b{Y,-e)-b{Y,_r\9)) 

y,H+l/2 


dr 


= I/O for alH < 0. Then, 
< ch{Ai + A2), 


where 
Ali = 


a 


-1 


^R+ 1/2 


dr 


and A 2 



a-^{b{Y,-,e)-b{yo-,e)) 

j.H+l/2 


dr 


Using the mean value theorem and the uniform boundedness of dxb, the second term is easily 
bounded as: 


A 2 < Ch 


Ys - yo\ 

gH- 1/2 ■ 


Similarly, we get that for all 0 < £ < 1, 


Ai < c sup 

0<r<s 



rpH—e 



ChS^ ^ sup 

0<r<s 



j^H—e 


Then, proceeding as in [29] (where Fernique’s theorem is invoked), we get the desired result. 

□ 


2.4. Contraction properties of the SDE. The inward Hypothesis 1.3 for the function b 
yields some contraction properties for the solution to equation (3), which will be useful for 
our purposes. We summarize some of them in this section. 

Let us start by giving a result quantifying the way two solutions starting from different 
initial conditions are getting close as t goes to 00 . 

Proposition 2.11. Assume Hypothesis (1-3). Let X be the solution to equation (3) start¬ 
ing from initial condition x, and Y the solution to the same equation starting from initial 
condition y. Then we have: 

\Xt-Y,\ < \x-y\e-^\ (25) 

where the constant a is introduced in Hypothesis 1.3. 

Proof. Recall that we have assumed; 

(6(x; 9) - b{y] 9), x - y) < -a\x - y\^. 

Now it is readily checked that the difference X — Y satishes 

Xt-Yt = {x-y)+ [ [5(X,; 9) - 6(W; 0)] ds, t > 0. 

Jo 
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Applying the change of variable formula in the Young setting (18), we thus get: 

\Xt-Yt\^ = \x-y\^ + 2 [\b{X,-,e)-b{Y,-,e),Xs-Ys) ds 

Jo 

< \x — y\'^ — 2a I \Xs — Ys\'^ds. 

Jo 

Our claim thus follows by a direct application of Gronwall’s lemma. 


□ 


The previous proposition can be enhanced in order to get a bound on increments for X — Y. 

Proposition 2.12. Under the same assumptions as in Proposition 2.11, the following in¬ 
equality holds true for 0 < u < v < oo: 

\5Xuv — dYuv\ < c\x — y\ — v\, (26) 

where the constant a is introduced in Hypothesis 1.3. 


Proof. The same computations as in the proof of Proposition 2.11 show that 

pv 

- 6Yu, = / [b{Xs; e) - b{Ys-, 9)] ds. 

J u 

We then use the mean value theorem and the uniform bound on dj), together with inequal¬ 
ity (25), which yields our result. □ 

We now study the contraction properties of the map h i-A- Y^ for h eH. 

Proposition 2.13. Assume Hypothesis 1.3. Let h be an element ofH and Y, Y^ he solutions 
to the system: 

yt = yo+ [ b{Ys]e)ds + ^ajBl, andY^^ = yo+f b{Y^^-, 6) ds+ ^ aj {Bl + hi) . (27) 
Jo Jo 

Then the difference Y^ — Y satisfies: 

< cexp ||h||^, (28) 

uniformly in t E M+. 


Proof. Consider s < t and proceed again as for the proof of Proposition 2.11 in order to get: 

\Y^h _ y^|2 < _2c, f -Yr\^dr + J2 f {^r " dj) dhi, 

Jo Jo 

and thus, applying Cauchy-Schwarz inequality we end up with: 


Wh _ W |2 

I s ■'si 


< -2a [ \Y^^ 

Jo 


d 

Yr\^dr + E \(Jj\ ||h-^'||vi/i,2([o,t]) 
i=i 



1/2 
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Now recall that ab < -^ + for a,b,e > 0. Applying this inequality to = 2a and 
resorting to Proposition 2.2, we obtain: 



The last inequality being valid uniformly in s < f, our claim is thus easily deduced by a 
direct application of Gronwall’s lemma plus the fact that exp(—a:f/2) is bounded over 

the time interval M_|_. □ 

We now push forward our analysis to include the increments of processes. 

Proposition 2.14. Assume Hypothesis 1.3. Let h he an element ofTL and Y he solutions 
to the system (27). Then for 0 < u < v, the difference Y^ — Y satisfies: 

- ^Yuv\ < Cl exp (-^) \\h\\n 

Proof. Like in the proof of Proposition 2.13, we start from the inequality: 

pv ^ pv 

\5Yi - 6Yu.f < -2a / {Y,^ -Y^fds + J^ {^s - U, <y,) dhf 

J u J u 

and since the first term on the right hand side above is negative, we get: 

d pv ^ 

\5Yi - ^ / {Y^^ - n, a,) dhi < sup in" - n| |a,|||h^||i_.ar,M]- 

j=i du se[M,q 

We now plug relations (10) and (28) in order to get: 

< Cl exp (-on) ||h||^(n - m)"". 

□ 

As a consequence of Propositions 2.13 and 2.14, the following bounds for the 7^-norm of 
the Malliavin derivative of the process and its increments hold. 

Corollary 2.15. Assume Hypothesis 1.3. Let Yt be the solution to the system (3). Then, 
for all t > 0, Yt belongs to and the Malliavin derivative satisfies that 

llA’nik < c exp , (29) 

uniformly in t ^ R_|_. Moreover, for 0 < n < n, 

\\D{5Y,,fi)\\ft < Cl exp {v - (30) 

uniformly in u and v. 

Proof. The fact that Yt belongs to follows from (28) and [34, Lemma 2, page 72]. In 
order to check (29), we first approximate Yt by a sequence of cylindrical random variables, 
use the definition of the Malliavin derivative for cylindrical random variables (see e.g. [34, 
Section 2.2]), and apply the inequality of Proposition 2.13. Then passing to the limit by 
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closability and taking "H norms, we obtain (29). Finally, (30) is proved similarly using 
Proposition 2.14. □ 

2.5. Ergodic Properties of the SDE. We recall here some basic facts about the limiting 
behavior of equation (3), mainly taken from [12]. We still work on the Wiener space {B, H, P) 
introduced in Section 2.1.1, and seen as the canonical probability space. Together with the 
shift operators 6^4 : 12 —)■ 12 defined by 

6tU!{-) = u{- + t)— uj{t), t G M, a; G 12, 

our probability space defines an ergodic metric dynamical system, see e.g. [14]. In particular, 
the measure P is invariant to the shift operators 6t, i.e. the shifted process (i?s(0r))sGR is 
still an m-dimensional fractional Brownian motion and for any integrable random variable 
F : 12 —)■ R we have 

lim - / F{et{u))dt = E[F], 

r^oo T Jq 

for P-almost all cu G 12. 

Under our coercive hypothesis on the drift coefficient of equation (3), the following limit 
theorem is borrowed from [12, Section 4]. Notice that its proof is based on contraction 
properties for the stochastic equation. 

Theorem 2.16. Let Hypothesis 1.3 hold, and consider the unique solution Y to equation (3) 
as given in Proposition 2.6. Then there exists a random variable Y : 12 —)■ R*^ such that 

lim \Yt{cu)-Y{etUj)\ = 0 

t—>-oo 

for P-almost all cj G 12. Moreover, we have E[|U|^] < cx) for all p> 1. 

It is worth observing that the convergence of Y towards Y can also be quantified in Holder 
norm, like in Section 2.4: 

Proposition 2.17. Let a G {0,H). There exists a random variable Z admitting moments 
of any order and a constant c > 0 such that for all 0 < s < t we have: 

\Yt-Yt\< Ze-'^^ and |(5 [U - F]< Ze"“(t - s)“. 

Proof. The proof is done along the same lines as for relation (26), and details are left to the 
patient reader for sake of conciseness. □ 

We now recall, as in [27], that the integrability of Y implies the ergodicity of equation (3): 

Proposition 2.18. Assume Hypothesis 1.3 holds true. Then for any 6 ^ G 0 and any f G 
C^(R'^;R) such that 

\fix)\ + \d^fix)\ <c{l + \xf) , xgR'^, 
for some c> 0, N ^ N, we have 

hm-/ f{Yt)dt = E[f(Y)], P-a.s. (31) 

T^OO T Jq 
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3. Proof of Theorem 1.7 


We can now gather the information we have obtained on the system (3) in order to complete 
the proof of Theorem 1.7. This will be divided in different steps. 


Step 1. The hrst step of the proof consists in applying Girsanov’s theorem. Namely, fix 
6 ^ G 0, and set 6 ^ = 0 + Then Proposition 2.10 implies that 


log 


(dVl 


I 


dPl 


\a-\[D^-^/%iY;er)]t - [D^-^/%{Y;eM^dt. 


1 

2 


(32) 


Step 2. We next linearize relation (32). To this aim, recall that b is dehned at Hypothesis 1.3 
as fe = deb. We add and snbstract the d-dimensional vector 

-ff) = 

where we abbreviate {[D^~^^‘^b(Y ; 9)]t, into [D^~^^‘^b{Y ; 9)]tU for notational sake. This 

easily yields: 


log 


dPg, 

dPl 


— h — I2 — -h — I4, 


(33) 


where 


Ii = ^l {a-^[D^-^^\Y;9)]tU,dWt)-:^ I \a-^[DX~^'"b{Y;9)]tu\‘^dt 


-irn^-i/2i 


'0 

I 2 = I {a-\[D^-^^h{Y- 9^)]t - [DX~^^'%{Y- 9)]t - [Dt-^''%Y- 9)]t{9r - 9)), dWt) 


2t 




H- 1 / 2 ; 


h = 


a-\[Dl-^'%{Y-9,)]t - [Dl-^/'%{Y-9)], - [Dl-^<'%{Y-9)]t{9, - 9))\^dt 


W-1/2 


H- 1 / 2 ; 


and 


Jo 

- [D^-^^%{Y- 9)],{9, - 9)), a-^[D^-^'%{Y- 9)\{9, - 9))dt. 


Step 3. In this step we set the ground for the identihcation of the main contribution to our 
log-likelihood. That is, we wish to show that as r —?■ cxd we have: 

u^T{9)u, (34) 

where r(0) is given by (6). Observe that (34) together with the multivariate central limit 
theorem for Brownian martingales (cf. [22, Theorem 1.21]), imply that Ji is the term that 
contributes to the limit in (5). In Steps 7 and 8 below we will show that I 2 ,l 3 and I 4 are 
negligible terms with respect to this main contribution. 
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In order to prove (34), according to the expression of we thus have to find an 

equivalent for the quantity 


Jr{Y)= / \a-^Nt{Y)\^dt, with Nt{Y) = 


ib{Yf,e)-b{Yt.r;e))u 

^H+l/2 


dr. 


(35) 


Furthermore, in order to dehne Nt{Y) properly, we have to extend the dehnition of Y to M_ 
by setting Yt = yo for all f < 0. Next we decompose Nt{Y) into Ni^t{Y) + N2^t{Y), where: 




{b{Yt;9)-b{Yt_r;0))u 

j.H+l/2 


dr, N 2 ,t{Y) = 


{b(Yt;9) - b{yo]9))u 

^H+l/2 


dr. (36) 


We shall bound those terms separately. 

Step 4- This step is devoted to reduce our computations to an evaluation for the expected 
value. Specifically, we shall prove that: 


lim -Eg [|W-E,[J.]|] = 0, 

T—>-oo 7" 


(37) 


where we write Jr instead of Jt{Y) in order to alleviate notations for the remainder of this 
step. Observe that owing to the Poincare type inequality (16), relation (37) can be reduced 
to show that: 


lim 

T—¥00 


= 0 , 


(38) 


which is what we proceed to do now. 
Indeed, by the chain rule we have: 


DJr = 2 {a-^Nt,D{a-^Nt))dt, 


where we have used the notation Nt for Nt{Y). Therefore, by Cauchy-Schwarz’ inequality, 
we obtain: 




cr 




'W| \\D{a-^Nt) Wndt. 

JO 

Thus, another use of Cauchy-Schwarz’ inequality yields: 


E<. [\\DJr\\n\ r 1 " 11 ^ II«IIL(M)_ 


1/2 




r r 

We shall prove (see the forthcoming inequality (46)) that E 0 [Jt-] is of order r. Hence: 


Y.g[\\D.Jr 


m < 


E 


1/2 


|cr (N.) \\-fl\\\2(^[0^r]) 


T 


1/2 


E 


1/2 


< 

r\j 


\Y> {N.) ||Hlli 2 ([ 0 ,.r]) 


- 1/2 


(39) 
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Let us now compute DNt- We can decompose W as in (36) and recall that Uq is considered 
as a constant, in order to get DNt = -DW.t + DN 2 ,t, with; 

rnkYuO), DYt) - {dMYt-r-,e),DYt.r))u 


DN,t = 


DN2,t = 


j-77+1/2 


dr 


(dMyuO), DYt)u 

^H+1/2 


dr. 


(40) 


We treat again those two terms separately. 

The term DN 2 ,t in (40) can be explicitly computed as {dxb{Yt; 9), DYt)u. Thus, 

according to Hypothesis 1.3(iii), we get: 


E 


1/2 


\D (A^2,-) Iklli2([0,r]) 


<E, 


1/2 


UO 


\DYtf^t^-^^dt 


By Corollary 2.15, is uniformly bounded in t > 0. This yields: 


E 


1/2 


\DiN 2 , 


l7dlL2([0,r]) 


A-H 


rV 2 

which tends to zero as r —?■ 00 since H > ^. 


< 


T 


1/2 


^H-l/2 ’ 


We next treat the term DNi^t in (40). Observe that it can be expressed as: 

/■* {{dMyf, 0), DYt - DYt-r) - {dMYt-r-, 0) - OMYu 0), DYt-r))u 

Jo ^ 77 + 1/2 

Therefore, by the mean value theorem and Hypothesis 1.3(iii), we get that 


dr. 


\\DN,Ah < 


\DY-DY 


t-r\\H 


j.H+1/2 


dr + 


\Yt r-Yt\\\DYt. 


r\\H 


j.H+l/2 


dr. 


We now appeal to Proposition 2.8 and Corollary 2.15, to obtain that for all e G (0, H), 
II \\D (W,-) lkllL2([0,T]) 


< 


dt 


exp (—c(t — r)) 
rW+i)/2 


dr 


+ zt 


dt (1 + 1 ) 


4e 


exp (—c(f — r)) 

/pl/ 2 -\-£ 


dr 


( 41 ) 


where is a random variable with finite moments of all orders. The integrals in the right 
hand side above can be bounded as follows: for f3 G (0,1) and a generic constant c > 0 which 
can change from line to line, observe that; 




_ /■* exp (—c(f — r)) 




dr = 
< 


exp (-c(t - r)) 
rP 


dr + 


exp (—c(f — r)) 


dr 


't/2 


exp{—cat)t^ ^ + t ^ Yt 

uniformly in t. This bound is useful for f > 1. When f < 1, we can bound by a constant c. 
Therefore, can be uniformly bounded by c(l + f)“^. Reporting this information into (41), 
we obtain that for r > 1 , 

rr 1 rr (i + t)4£ 


l|-D^l,-|l'HllL2([0,r]) 


< 


(1 +^ 7o (1 + ^) 


dt < 


l+2e ~ £ 
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Thus, for r > 1, 


E 


1/2 




hNi^Uo.t]) 


- 1/2 




which tends to zero as r —?■ cxd as long as we choose £ < Gathering our bounds on DNi, 
DN 2 and plugging them into (39), the proof of (38) is now completed. As a consequence of 
this inequality together with the Poincare type bound (16), we get that: 

Jt — Ee [Jr 


lim Efl 


= 0 . 


This allows us to be reduced to an evaluation of Eg [Jr] only in the sequel. 

Step 5: We now show that one can approximate the quantities related to Y by quantities 
related to the stationary solution Y . Namely, we denote by Jr{Y) the same quantity as Jr{Y) 
(see (35)), with Y replaced by Y, and we set Nt(Y), N2^t(X) accordingly. Observe 

that, with (34) in mind, we have to compare Nt{Y) with Nt{Y) up to some 0(1) terms. We 
will achieve this aim thanks to our decomposition (36). 

Step 5a. The term A" 2 ,t(E) is easily estimated; it is readily checked that 

cniK^uO) - b{yo;9))u 


N 2 AY) = 


tH- 1/2 


Invoking Proposition 2.7 as well as the linear growth for b, we get hmt_,.oo N 2 AY) = 0, almost 
surely and in any U’iyi) for p > 1 . 


Step 5b. Let us compare Ni^Y) and W,i(E). One can write: 


n.ay) - NiAY) = 

where we have set: 




t—r,t 


u 


j.H+l/2 


dr = Wi,t + Ni 2 ,t, 


Niu = 


^t/2 5[b{Y-e)-b{Y-e)\,_^^^u 


j.H+ 1/2 


dr, 


Nl2t — 




It/2 


j.H+1/2 


dr. 


We now show that those two terms give raise to a negligible (that is o(l) in terms of t) 
contribution. 

The term involves some increments oi Y — Y, and should thus be handled thanks 
to Proposition 2.17. However, due to the composition with the nonlinearity b, we have to 
resort to a special treatment concerning rectangular increments. To this aim, let us define a 
path a : [ 0 , 1 ]^ —?■ M of the following form; 

d) — Yt-r + A(17 — Yt-r) + y{Yt-r ~ Yt_r) + A/i AH 

where we have set ^Yt-r,t = Yt — Yt-r — Yt + Yt-r- We also write A for the composition of 
a with the non linearity b, that is A{X, p) = b{a{X,fi)). It is then readily checked that: 

b{Yt)-b{Yt-r)-b{Yt) + b{Yt-r) = A(1,1 )-A(1,0)-A(0,1) + A(0,0) 

= / d\^A[X, p) dXdy. 


(42) 
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Furthermore, thanks to the fact that the hrst and second order derivative of b are bounded 
(see Hypothesis 1.3), we get that: 

fi)\ < \dxa{X, fi) df,a{X, fi)\ + \dl^a{X, fi)\ 

< + I {\Yt - Y\ + I + I AF,_,,i| 

< (|<5F-.,t| + IAF_,,i|) (|Fi - Ftj + \Yt_r - Yt-r\) + \AYt-r,t\ 

We now bound 6 Yt_r^t thanks to Proposition 2.8 and Yt — Yt, Yt-r — Yt-r and AYt-r,t 

according to Proposition 2.17. Gathering all those estimates we get that for all a < and 

e < H, the following holds true: 

\dl^A{X, /i) I < (Z (1 + + Z Z + Ze-"(*-"V“. (43) 

Therefore, reporting (43) in (42) and then in the definition of and taking a > H — ^ 

and e < i, we get that: 


Nii,t < Z 


fb2 ^-c(t-r) 


1‘1 — CL 


— dr + Z(1 + f) 


2e 


'•p2 ^-c{t-r) 


n/ 2 +e 


dr 


< + z{l + 

It is now readily checked that limt_,.oo iVn ^ = 0 almost surely with uniform L^{Q) bounds, 
and thus yields a negligible contribution. 

The term iVi 2 ,t is bounded in the same way for f > 1: choosing a > H — 1/2 and £ < |, 
we have: 

ft ^-c{t-r) 


Nu,t < ^ 


rt -c{t-r) 

— dr + Z{l + tr 


< 

r\j 


/i/2 

z 


/ 2 — o. 


h/2 


fy\ j2-\-S 


dr 


e-d-r) t dr < 


j-H+l/2-a tl/2+e 

which gives another negligible contribution. 


+ 


/i/2 


^H+l/2-a ' 


Step 5c. We should also bound the term Notice that Y is not constant on M_, so 

that: _ _ 

We shall simply bound this quantity as: 


l^ 2 ,i(H)| < \u 



\h{Yt-e)\ + \h{Y,_r.e)\ 

^H+l/2 


dr. 


(44) 


In addition. Proposition 2.8 can be shown for Y exactly in the same way as for Y. Therefore, 
for any e > 0 we get: 

UYp e)\ + \b{Yt_r; d)| < Z, (1 + + |r|2-), 

where Z^ is a random variable admitting moments of all order. Plugging this bound into (44), 
we get a bound of the form |iV 2 ,i(y)| < for an arbitrarily small constant e. 

Summarizing Steps 5a to 5c, we have found that Nt{Y) — Nt{Y) is of order t~'^Z with 
r] > 0 and Z G np>iLP(r2). Hence Jt{Y) — Jt{Y) is of order which is a negligible term 
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on the scale r. This allows us to consider the limiting behavior of Jt{X) instead of JriY)^ 
which is what we do in the sequel of the proof. 


Step 6. We are now reduced to the analysis of the quantity E 0 [JT-(y)], which is equal to 
u^'^u where: 


T = 


dt 


Ee[{biYp,e)-b{Yt-rP,e)y 


UT 


-DT^-l 


yb{Yp,9)-biY,.r,;e))] 




r 


l/2+H 1/2+H 


dridr2- (45) 


To 


By stationarity of T, the expected value does not depend on t\ 

E4(6(F,i 0) - ; »))T(ff-')T(r-i(i(F,; 0) - »))] 

= E„[(6(F„i 0) - KF„i «))’'(ff-i)Tff-i(j(F„i 0) - b{Y,„ «))|. 


Furthermore, we have already observed that Proposition 2.8 is also satished by Y. Hence, 
thanks to the fact that 6 is a Lipschitz function in x, the following bound holds true: 

|E„[(S(F„i«) - S(F,.i«))T(ff-i)Ta-i(S(F„; 0) - i(F,,; »))|| < (r« A 1) (r" A l) , 

Plugging this information into (45), we obtain that 4/ is a convergent integral and: 

^e[Jr(Y)] = ™^r(0)M, (46) 

where r( 6 ') is the convergent integral given by ( 6 ). In particular, we have: 


lim 

r—>-co 7 " 


u^T{e)u. 


Summarizing our considerations up to now, we have proved (34). Hence Ji in (33) converges 
to the term m'^A/'(0, r(6*)) — ^u^T{9)u in (5). 

Step 7. We next show that the term J 3 in (33) converges to zero in Pe-probability as r ^ cxo. 
For this, set gt{9) = a~^[Dy ^^‘^b{Y;9)]t. One can recast J 3 into: 


h = r Igtidr) - gtiO) - dggy9){9r - 9)^ dt, 

Jo 

and we also recall that 9^—9 = In addition, a simple application of Taylor’s expansion 

for multivariate function yields the existence of a A G (0,1) such that: 


gt{dr)-gt{Q) = dQgt{i>){9r-9), where ix = d Y \{9r - 9). 


Notice that under our standing Hypothesis 1.3 we have dggt{^x) = a ^^‘^b{Y ; ^a)]*- We 

thus get: 


h = — W ^MtiY^u^dt, where MfCY) = 

^ Jo 


5[b{Y;^x)-b{y;' 




J.I/2+H 


dr. 


Next we decompose MtiY) into Mi^tiY) + M 2 ^t{Y), where 


M,,t{Y) = 


6[b{Y-,^x)-b{Y-,9)]t-r,t 

J.1/2+H 


dr, M 2 ,t{Y) = 


5[b{Y-,ix)-h{y-.0)]t-r,t 

J.1/2+H 


dr. 
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Recall that we have extended the definition of Y to M_ by setting Yt = i/q for all f < 0. 
Thus, using Hypothesis 1.3, we obtain that 

T-' r \a-^M2t(Y)u\-^dt<T-^- 

Jo Jo ^ 

Thus, by Proposition 2.7, the L^(r2)-norni of this term is bounded by . Hence, this 

term converges in P^-probability to zero as r —)■ cxo. 

On the other hand, fix a G (’ 271)’ write 


M.AY) = 


(< 5 [fe(y; 6 ) - Hy; e)]t-r,tnmy-, a) - Hy-, o)]t-r,ty-‘^ 


J.I/2+H 


■ dr. 


Then, appealing to Hypothesis 1.3, we get that 


mAy)\ 


< 


1 ni + \Yt\'^ + \Yt_rn\Yt-Yt_, 


I 1 —Q 


T 


■a/2 


[1—a) j^l/2-\-aH 


■ dr 


Therefore, 

\M,Ay)\" 


t rt 


< 

r\j 



0 JO 


(1 + + \Y^_rT)\Yt - (1 + 1 ^* 1 “ + \Y^_rX) % - yt- 

H{l-a) 1/2+aH H{l-a) l/2+aH 

1^1^ Inin 


1 1—a 


r2 I 


dridr2. 

1 ' 1 ' 2 ' 2 
Thus, again by Proposition 2.7, the L^(r2)-norm of the term - J/ \a~^Mi t{Y)u\^dt is bounded 
by 
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7-1+a 


T ft Pt 



0 Jo Jo r 


l/2+aH 1/2+aH 


dridr2dt = 


/ ,l-2Ha 


-l+a 




Ch 


^a{2H+l)-l ’ 


which converges to zero as r —)■ oo since a > 


2H+1 ■ 


Step 8 . We finally show that I 2 and 14 are also negligible terms. Since is the quadratic 
variation of the martingale I 2 , this implies that I 2 converges to zero in P^-probability as 
r —)■ 00 . On the other hand, applying Cauchy-Schwarz inequality to J4, we get that 

I/ 4 I < \a-^[D/-y\Y-e)]tu/d^ ' x/ 3 '/^ 

Thus, by the results in Steps 3 and 4, we obtain that J 4 converges to zero in P^-probability 
as r —?■ CX 2 , which completes the proof of the theorem. 
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